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PREFACE 


The subject “Recursion Sequences” is related with the secon¬ 
dary school curriculum (arithmetic and geometric progres¬ 
sions, the sequence of squares of natural numbers, the sequence 
of coefficients in the quotient of two polynomials in 
ascending powers etc.). It is at the same time a real 
though a small mathematical theory*, complete and clear 
like everything created by the eminent masters of analysis 
who founded it. 

The fundamentals of the theory of recursion sequences 
were worked out and published in the twenties of the 
XVIII-th century by the French mathematician Abraham 
De Moivre [the formula (cosa-f i sin a) n = cosna + i sin na 
bears his name] and by one of the first members of the 
Petersburg Academy of Sciences, the Swiss mathematician 
Daniel Bernulli. A comprehensive theory of these sequen¬ 
ces is due to the great XVIII-th century mathematician, 
member of the Petersburg Academy academician Leonard 
Euler who devoted to the recursion sequences (series) 
Chapter 13 of his Introduction to the Analysis of Infini¬ 
tesimals (1748). Among the recent works the courses of 
the calculus of finite differences read by the renowned 
Russian mathematicians academicians P.L. Chebyshev and 
A. A. Markov should be especially noted; the courses 
include notes on the theory of recursion sequences. 


* For the reader who is versed in analysis it can be stated that 
the theory is a strict analogue of the theory of linear differential 
equations with constant coefficients. 
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1. The concept of a recursion sequence is a broad gene¬ 
ralization of the concept of an arithmetic or geometric 
progression. Its particular cases also include the sequences 
of squares or cubes of natural numbers, the sequences of 
digits of a rational repeating decimal (and any recurrent 
sequence in general), the sequences of coefficients in a 
quotient of two polynomials in ascending powers of %etc. 
It follows that in the course of mathematical studies at 
a secondary school recursion sequences are very often met 
with. The theory of recursion sequences is a special part 
of the so called calculus of finite differences. This theory 
will be expounded here for readers with no previous 
experience in the subject (in only one place we shall 
refer without proof to a general proposition established 
by the theory of linear algebraic equations). 

2. Sequences will be written in the form 

Uj, u 3 , u 3 , ..., u n , . . . , (1) 

or, shortly, {«„}. If there is a natural number k, and 
numbers a u a 2 , ..., a k (real or imaginary), such that 
starting from a certain number n and for all following 
numbers 

+ a 2 «„ + A_ 2 + • • • +a k u n (n>m>l), (2) 

the sequence (1) is called a recursion sequence of order k 
and the relation (2), a recursion relation of order k. 

Thus, a recursion sequence is characterized by each of 
its terms (starting from a certain one) being expressed by 
the same number k of terms immediately preceding it, 
according to formula (2). The very name “recursion se¬ 
quence” is used just because one must turn back to the pre¬ 
ceding terms in order to obtain the following one. Some 
examples of recursion sequences will be given below. 

Example 1. Geometric Progression. For the geometric 
progression 

«i = c, « 2 = a<7» m s = «<7 2 > •••. u n = aq n ~\ ...; (3) 

relation (2) takes the form 

U n + l = QUn ( 4 ) 

where k — 1 and = q. Thus the geometric progression is 
a recursion sequence of the first order. 
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Example 2. Arithmetic Progression. For the arithmetic 
progression 

u 2 = a-\-d, u 3 = a-\-2d, ..., w„ = a + (n — \)d, ... 
we have 

u n +1 = u n A~d. 

This relation is not of the form of relation (2). Characte¬ 
ristic of the latter is that its right-hand member contains 
only sequence terms with constant coefficients. However, 
considering two relations for two consecutive values of n: 

“». + * = “« + i + d and = + d 

and subtracting the latter from the former term by term 
we obtain 

U n + 2 u n + l ~ U n + ] 

or 

U n + 2~ 2« n + l U rr ( d ) 

Relation (5) is of the form of (2). This suggests that the 
arithmetic progression is a recursion sequence of the second 
order. 

Example 3. Consider the ancient problem of Fibonacci * 
about the number of rabbits. Required is the number of 
pairs of adult rabbits descending from one pair in the 
course of a year if it is known that each adult pair of 
rabbits gives birth to a new pair every month, the new¬ 
born becoming fully adult in a month. The interest of 
the problem lies not in its solution, which is easily arrived 
at, but in the sequence whose terms express the total 
number of adult pairs of rabbits at the initial moment 
(«i), in a month (u 2 ), in two months ( u 3 ) and in general 
after n months (u„ + 1 ). Obviously u 2 — 1. One newborn pair 
will be added a month later, but the number of adult 
pairs will be the same: « 2 =1. In two months (from the 
initial moment) the little rabbits will become adult and 
the total number of adult pairs will equal 2: u 3 = 2. Assume 
that we have calculated the number of adult pairs after 


* Fibonacci, or Leonardo da Pisa, a medieval Italian mathematician 
(circa 1170), wrote the book Liber Abaci, containing extensive arith¬ 
metic and geometric information borrowed from peoples of Middle 
Asia and Byzantium and creatively reworked and developed by him 


8 




n —1 months, u,„ and after n months, «, 1 + 1 . By this time 
the u n adult pairs that existed before will have an off¬ 
spring of u„ additional pairs, so that after n-|-l months 
the total number of adult pairs will be 


Hence 


^n + 2 ^n + 1 U n . 


( 6 ) 


w 4 = u 3 + « 2 —3, = u 4 -f « 3 == 5, 

« 0 = « 5 -f m 4 ==8, « 7 == «„-I- « 8 = 13, .., 


Thus the sequence obtained is 

u 4 = 1, u 2 = 1, w 3 = 2, « 4 = 3, 1 
w 5 — 5, u 6 — 8, — 13, ••*, J 


i.e. every following term is equal to the sum of two pre¬ 
ceding terms. This sequence is known as the Fibonacci 
sequence (series) and its terms, as Fibonacci numbers. 
Relation (6) shows that the Fibonacci sequence is a recur¬ 
sion sequence of the second order. 

Example 4. Consider now the sequence of squares of 
natural numbers: 

l 2 , « 2 = 2 2 , u 3 = 3 2 , ...,«„ = n 2 , ... (8) 

where u n+1 = (n 4-1) 2 = n 2 + 2/1 + 1. Thus 

«» + i = “« + 2 a + l. (9) 

Increasing n by unity we have 

u «+! = u n+i + 2ft-f3. (10) 

And hence [subtracting (9) from (10) term by term] 

+ 2 ~~ +1 = U n + i +1 + 2 , 

or 

U n + 2 = 2« n + i U n - (-2. (11) 

Increasing in equality (11) n by unity gives 

u «+3 = 2 m„ +2 w„ +1 + 2. (12) 


Hence [subtracting (11) from (12) term by term] 
u n+3 u n+2 — 2u n+2 3 u n+1 + u n , 


or 


( 13 ) 


2-257 


u n + 3 — 3 u n + 2 3«„ +1 + «„. 
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The recursion relation obtained is of the third order. 
Hence sequence (8) is a recursion sequence of the third 
order. Similarly it can be verified that the sequence of 
the cubes of natural numbers 

I s , 2 :! , 3 3 , .... n\ ... (14) 

is a recursion sequence of the fourth order. Its terms 
satisfy the equation 

u n+i ~ 4«„ +3 6u„ + 2 -|- 4 u n+1 u n , (15) 

and the reader is invited to derive it. 

Example 5. All recurrent sequences are recursion sequences. 
Consider for instance the sequence of digits of the rational 
repeating decimal 

^ = 0,57132132132 ... 


In this sequence 

Uj = 5, u 2 = 7, u 3 = 1, « 4 = 3, ) (16) 

u 6 = 2, u 6 = 1, u 7 = 3, ... I 

It is obvious that 

= (n> 3). (17) 

Rewriting equation (17) to give it the form of relation 
(2) we have 

u n+3 = 0-u n+2 -\-0-u n+1 -f- 1 • u n . 

Thus, Eq. (17) is a recursion relation of the third order 
(k — 3, a x =0, a 2 = 0, a 3 =l), and sequence (16) is a recur¬ 
sion sequence of the third order. 

Example 6. Now let us consider the sequence whose 
terms are the coefficients of the quotient of two polyno¬ 
mials arranged in ascending powers of x. Given 

P (x) — A o -f~ A x x + ■ • • + A L x l 
and 

Q (x) = B 0 - j- B x x + . .. -f B k x k (B 0 =4= 0). 

Let us divide P (x) by Q (x); if there is a remainder 
the division can be infinitely continued. The quotient will 
contain terms in the following order: 

D 0 + D x x -j- D a x a -j- D a x 3 + ... + D n x n + • • • 
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Consider the sequence 

^ 1 = D 0 , u 2 — D^ } . .., u n = D n _i, ... (IS) 

We shall prove that this is a sequence of order k (remem¬ 
ber that k is the power of the divisor). In order to do it 
we choose an arbitrary natural number n satisfying only 
one condition —fe-j-1, and in dividing we stop at 
the term of the quotient with x n+k . The remainder will 
be a polynomial R (x) containing powers of x higher than 
n + k. Writing down the relation between dividend, divisor, 
quotient and remainder we obtain the following identity: 
A 0 -f • • • ~\~ AiX l = 

= {B 0 J r ■ ■ • + B k x k ) • (D 0 + . .. + D n+k x n+k )-\- R (x). 

Find the coefficients of x n+k in the right-hand and left- 
hand members of the identity and equate them. As 
n + k^l + 1, the coefficient of x n+k in the left-hand member 
is zero. Hence the coefficient of x n+k in the right-hand 
member is also zero. But the terms with x n+k are present 
only in the product (£„+ . . . + B k x k ) • (D 0 +... +D n+k x n+k ). 
(The remainder R (x), as was mentioned above, contains 
only higher powers of x). Thus the required coefficient is 

D n +hB 0 -i-D n + k _ 1 B 1 -j-D n B k . (19) 

According to the foregoing argument this coefficient is 
equal to zero: 

D n +kB o + B) n+k ^ 1 B 1 + • • • -f D n B k = 0. 

Hence (as B 0 ^=0) 

= (n>l-k + 1). (20) 

This is a recursion relation of order k. It follows that 
sequence (18) is a recursion sequence of order k. 

3. Of all the above examples Example 6 is of the most 
general character. We shall now show that any arbitrary 
recursion sequence of the k-th order 

Ult U^y * * • , U n , • • > , (21) 

satisfying the equation of the form 

“n+k = aiU n+ k-i+ ■■■+a k u n (n>m>l), (22) 


2 * 
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coincides with the sequence whose terms are the coeffi¬ 
cients of the quotient of a polynomial P (x) by the poly¬ 
nomial 

Q (x) == 1 — a x x — . .. — a k x k . (23) 

Let n be an arbitrary natural number satisfying the 
condition n>ft + m — 2; multiplying the polynomial Q (x) 
by Uj-f u 2 x + u 3 x 2 -j- + X x n we obtain 

(1— a x x—a 2 x 2 — . . .—a k x k )(u 1 + u 2 x + 

+ • • • + u li + m-l xk + + • • • + “n +1*”) — 

= [ U 1 + ( u 2 — a i u i) x -p • • • + 

-f- ... *f- (u k + m -i 1 - d x U k + m _2 • • • d k U m _ x ) X k + m 2 ] -f- 

+ [(« k +m — a 1 u k+m _ l —...—a k u,„)x k+m - 1 ) + 

+ • • • + (Un + i—a 1 u n —...—a k u n _ k+1 )x n ] — 

— [( a i u «4l + • • • + a A u n-ft + 2 )*' I + 1 + ■ ■ ■ -\- a k u n + iX n + k ]. (24) 

The polynomial in the first pair of square brackets is 
of a degree not higher than l — k + m —2; the coefficients 
in it are independent of the chosen number n. Let us 
call it P(x): 

P(x) = u 1 + (u i —a l u l )x -f 

+ • • • + (^A + m-l 4/71-2 • • • a k U m - 1) X k+m 2 . (25) 

All the coefficients in the polynomial in the next pair 
of brackets are zero by Eq. (22). Lastly, the third pair 
of brackets embraces a polynomial whose coefficients de¬ 
pend on n\ it contains no terms of a power lower than 
rt-f 1. 

Denoting this polynomial by R n (x) we can rewrite 
identity (24) in the following form: 

P (x) = (1 — a x x — a 2 x 2 —... — a k x h ) (u x -\-u 2 x-\- 

+ • • • + u n + l x ") + Rn ( x )- (26) 

Hence, it is obvious that u t + UgX+ ... +u n+1 x 11 is the 
quotient and R n (x), the remainder of the division of P (x) by 

Q ( x) — 1— a x x — a t x 2 — ... — a k x k 

i. e. 

dy, w 2 , •. •, u n , u n + j, ... 

is indeed the sequence of coefficients of the quotient of 
polynomial (25) by (23). 
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Consider, for example, the Fibonacci sequence 
= 1, u 2 — 1, u 2 = 2, u £ - o, 5, . *. 

Since its terms satisfy the relation 

“n + 2 = tJ n+1 + u„ (n> 1), 

then m=l, /a = 2, = 1, n 2 =l and Q (x) = 1— x—x 2 . 

The highest power in P(x) must not exceed k-\- m — 
— 2=1. By formula (25) we have 

/>(*)=!+(! — 1 x 1) jc= 1. 


Thus, the Fibonacci numbers coincide with the terms 
of the sequence of coefficients of the quotient of 1 by 
1 — x —x 2 . 


4. One of the problems which are offered for solution 
to the students of secondary schools when studying the 
arithmetic and geometric progressions and the sequence of 
squares of natural numbers is to find the sum of n terms 
of each of these sequences. 

Let 


u lt ll 2 , 


u 


n> • • • 


(27) 


be a recursion sequence of order k whose terms satisfy 
the relation 

Un +k = a 1 u n+ri _ 1 + a 2 u tl+k _ i +...+a k u n (n^m). (28) 

Consider a new sequence formed by the sums s„ of 
numbers (27); 

$, = «!, s 2 = «! + «„ ..., s B = u 1 + u l + ... +u„, ... (29) 

We shall show that the sequence of these sums is also 
a recursion sequence of order k +1 and its terms satisfy 
the relation 


S «+A+l — (1 + fl l) S « + ft + 

+ (<*.—«j) *»+*-!+ • •• +(«*— a k s n . (30) 


Note that 


Ui —Sj, — S 2 Hi — ^2 ^i» ••• 

• • •. «»=■»„—(«! + • • • +« n -i) = S,— S„-i, 
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Assuming s o = 0 so that u 1 = s 1 — s 0 and substituting in 
equation (28) for u x , u 2 , . .., u n their expressions in terms 
of s 0 , s lt ..., s n we obtain 

+ k S n + k-l (^n + k — l + It- 2 ) T~ 

"F ^2 (^n +A- 2 ^n + k- 3) “H • ■ ■ ~f~ (S„ 


whence 

S n + *= 0 + fl l) S h+A -1 + t a 2 fl l) S n + k -2 +••• + 

+ («*—o*-i)s n — a A s„_, (n>m), 
or,replacing n by n+1 we obtain 

s n + £ + l = 0 + a i) S n + * + ( a 2 a i) s « + A-l + ••• + 

+ ( a k—a k - 1 )s n+1 —a k s n (n>m —1). 

This is a recursion relation of order k-\-\. Several 
examples are given below. 

(a) Geometric Progression. In a geometric progression 

u n = aq n ~ 1 and s n = u x + u 2 + ... +u„ = a + aq + .. .aq n ~ x . 
Since the terms of { u„\ satisfy the relation u, l + 1 = qu n , 
the terms of {s„} must satisfy the relation 

«»+* = (! +9)s« +1 — qs a . (32) 

(b) The Sequence of Squares of Natural Numbers. Here 
u n — n 2 and s n — 1 +2 2 + • • • +n 2 . Since the terms of {«„} 
satisfy the relation 

W n + a = 3w n + 2 3m„ + 1 -h 

(cf. p. 9) 

s n+4 ~ 4 s n + 3 ^ s n + 2 “1“ 4s n + 1 

(c) Fibonacci Numbers. Since the Fibonacci numbers satisfy 
the relation 

U n + 2' ^n + i “F U n 

their sums s n must satisfy the relation 

s n + 3 ~ + 2 S n . 

5. In the case of the simplest recursion sequences, for 
instance, arithmetic and geometric progressions, sequences 
of cubes or squares of natural numbers, and recursion se- 
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quences, any term can be found without calculating the 
preceding terms. In the case of the Fibonacci number se¬ 
quence or of the general sequence of coefficients in the 
quotient of two polynomials it seems, at first sight, to 
be impossible to do that. In order to calculate, say, the 
thirteenth Fibonacci number we preliminarily find all 
the preceding terms (using the relation M„ +2 = M n+1 + u n ): 

U-I — 1 * U 2 — 1 * m 3 = 2, m 4 = 3, m 5 = 5, u 6 = 8, 
u 7 =13, m 8 = 21, m 9 = 34, m 10 = 55, m u = 89, 

M n =144, m 13 = 233. 

We now turn to a detailed study of the structure of 
recursion sequence terms. This study results in formulas 
which make it possible to calculate in the most general 
case any terms of a recursion sequence without resorting 
to the calculation of the preceding terms. These formulas 
can be regarded as far reaching generalizations of the for¬ 
mulas expressing the general term of an arithmetic or geo¬ 
metric progression. 

Consider the following recursion relation of order k : 

^n + k ~ + A-1 “t~ + "F • • • "F • (33) 

If it holds for any value of natural numbers n=l, 2, 
3, .... then assuming n=l we can obtain 

U A + 1— a i U ft + a 2 W *-l+ • • • +%«!• 

Hence, knowing u u m 2 , .... u k , we can now find u k+1 . 
Assuming n = 2 in relation (33) we obtain 

U A + 2 = fl l U 7; + X + a 2 W A+ • • • +<W 

So now we know the value of u k+ „_ too. In general if m 
is a natural number and we have calculated the terms 

M x , M 2 , ..., U/ t , llfy 4-1 , •••> U m + !l _ i 

of a sequence, then assuming n — m in relation (33) we 
find the next term u m+k . 

Thus the terms of a recursion sequence of order k which 
satisfy relation (33) are unambiguously determined by this 
equation if the first k terms of the sequence u x , u 2 , ..., u k 
are known. Choosing the sequences in various ways (there 
are no conditions limiting this choice) we can obtain an 
infinite number of different sequences which satisfy rela- 
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tion (33). The distinction between them will manifest 
itself in the first k terms (at least in one of them) and 
will be observed in the following terms too. 

For example, the relation of the first order 

u n +1 = c ! u n 

is satisfied by any of the geometric progressions whose 
common ratio is q (they differ in the value of the first 
term «*); the relation of the second order 

w n+2 ~ 2u„ +1 u n 

(or u n+2 — u n + l — u n+1 — u n ) is satisfied by any of the arith¬ 
metic progressions which differ in at least one of the terms 
u l = a and u 2 = a-\-d, i. e. differ either in the value of 
the first term (a) or in the value of the common diffe¬ 
rence ( d ) or in both. 

Now consider another relation of the second order 

^« + 2 = +1 T 

Besides the Fibonacci sequence 

1, 1, 2, 3, 5, 8, 13, 21, 34, ... 

where u 1 = « 2 =l this equation is satisfied by many other 
sequences obtained by sampling various values of u x and u 2 . 
For example, with u x = —3 and u 2 = 1 we obtain the fol¬ 
lowing sequence: 

—3, 1, —2, —1, —3, —4, —7, —11, —18, —29, ... 

Let us have a certain number of sequences which sa¬ 
tisfy the same relation (33) 


U\ i i Dn> • ■ • > 


Zlt ^2> • • • > Z n , * • * i 

The following equations will hold : 

X n + rt ==a i X n + r<-l+ a 2 X n+ft-z J r ■ ■ ■ +d/t x n> 
Un + k = ^li/n + k-l ~f~ ^iUn + k-t ~F • • • 


2-n + k ^l^n + k-1 ^2^n+/i-2 I • • • 
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Taking as many arbitrary numbers A, B, ..., C, as 
sequences (34), we multiply each term of the first of equa¬ 
tions (35) by A, of the second, by B, . .., of the last, 
by C and add the products: 

Ax,,*. s+ By n+k + ■ ■ ■ +Cz n+k — 

— a 1 {Ax n+h _ 1 -f By n+k _ l - J r ■ • ■ + Cz n+k _ 1 ) 4 
+ (Ax n+k _ 2 + By n+k _ 2 + ... -\-Cz n+k _ 2 ) -f- 
+ ■ • • J r a k {Ax n + By,, + • • • + Cz n ), (36) 


It follows that the sequence 

t 1 = Ax x -f- By x -f • • • + Cz lt 
t 2 = Ax 2 4" By a 4- • • • + Cz 2 , 


t n — Ax n 4- By n 4- • • • 4- Cz n , 


(37) 


obtained from sequences (34) by multiplying each term of 
the first sequence by A, of the second by B, ..., of the 
last by C and then adding the products termwise (first 
terms with the first ones, second terms with the second 
ones etc.) satisfies relation (33). Since the choice of the 
numbers A, B ..., C is arbitrary we can by changing 
them obtain, generally speaking, different values of the 
terms t lt t 2 , t s , ... . 

Now let the sequence 

«i, «*. •••,«„, • • • (38) 

be such that it satisfies relation (33); the question is: can 
we give the numbers A, B, . .., C the values such that 
the first k terms of sequence (37) coincide with the first k 
terms of sequence (38)? If this can be achieved then accord¬ 
ing to the preceding argument all the terms of sequences 
(37) and (38) will also coincide, i. e. we will have for 
any natural number n: 

u n = Ax„ 4- By„ 4- • • • 4 Cz n (39) 

Thus it becomes possible (as yet hypothetically) to ex¬ 
press any of the infinite number of sequences which satisfy 
the same recursion relation of order k in terms of some of 
sequences (34), by formula (39). The realization of this 
possibility depends on whether we can select the numbers 


3-257 
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A, B, ..C so that the equations 

Ax x -f- By 1 + • • • + Cz 1 — u 1 , 

Ax 2 + By 2 + • ■ • 4* Cz 2 — u 2 , 

Ax k + By k -f- . . . + Cz k = u k 

with arbitrary given values of their right-hand members u lt 
u 2 , be satisfied. 

Since the numbers A, B, . . ., C are in this case the un¬ 
knowns and the number of equations is equal to the £-th 
order of the recursion relation, it follows, that the number 
of the unknowns A, B, ...,C [which coincides with the 
number of sequences (34)] should also be taken equal to k. 
It is known that a system of k algebraic equations (40) 
in k unknowns A, B, . . ., C can have solutions depending 
on the kind of the coefficients in the system: x lt y lt ... 

. . ., z lt .. ., x k , y k , . . ., z k , i. e. on the kind of the first 
terms of sequences (34). The solutions will certainly exist, 
with arbitrary right-hand members u lt u 2 , .... u l{ if we 
assume, for instance, that 

x t = l, y x = 0, ..., z t = 0; 
x 2 0, y 2 = 1, - ■ ., z 2 = 0, 

*a = o. y* = o, ■■■, z k = \. 

Indeed, in this case system (40) takes a very simple 
form, and the solution is directly obvious: 

A =«!, 

B = « 2 , 

C = u A 

Of course a different choice of numbers is possible with 
which system (40) has a solution whatever the right-hand 
members of the equations. Assume, for instance, 

*i= 1, yi = 1. • • *1 

* 2 = 0, y 2 = 1, ..., z* 

*a = 0, y k = 0, .... z fc 
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The system now takes the form 

A -f- B + • • • + C — u lt 'j 

B + . • ■ + t = U 2t i 
• • • .. 

C = u k . 

Hence we obtain 

C — u r ., B — u 2 — u 3 , A=u t — u 2 ■ 

We now turn to the general case and state the following 
theorem: 

For the system of k linear algebraic equations (40) in k 
unknowns to have one and only one solution A, B, .... C 
with any values of the right-hand members u ,, u 2 , . . u u 
it is necessary and sufficient that the corresponding homo¬ 
geneous system 

Ax ! By 1 + • • • + Cz x = 0, 

Ax 2 -j- By 2 + • • • -f- Cz 2 = 0, 


Ax h -f By k -T • • • + Cz k — 0 

have one only zero solution : * 

A = B= ... =C = 0 

The reader will easily ascertain that the condition in 
the formulation of the theorem is fulfilled in the particular 


* This is a convenient proposition since its application does not 
require any knowledge of the theory of determinants. The reader who 
is familiar with the theory will remember that for system (40) to 
have a solution with any values of the right-hand members of these 
equations it is necessary and sufficient that the determinant of the 
system 



*i Vi ■ 

■ Zl 

A = 


■ 22 


x n y n ■ 

• z„ 


be different from zero. The same condition is necessary and sufficient 
for system (40) to have one and only one solution for any given 
right-hand members (say, equal to zero). Thus for a system of ft linear 
equations in ft unknowns the condition for the existence of a solution 
with any right-hand members coincides with the condition that one 
and only one solution with zero right-hand members exists. Just this 
fact is expressed by the statement in the text. 


3* 
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cases (41) and (42). Later on we shall meet with cases in 
which the proposition stated above will prove useful. Until 
then we shall simply make use of the fact (established 
irrespective of the theorem) that there alwavs exist numbers 
x lt . . ., z x , .... x k , . . ., z u [the first terms of sequence (34)] 
such that the system of equations (40) has a solution with 
any values of u lf « 2 , . . ., u k . 

If such numbers have been chosen as the first terms of 
sequences (34), then according to the foregoing any sequence 
which satisfies the recursion relation (33) can be expressed 
by formula (39) in which the numbers A, B, C are 
found from equations (40). The system of k sequences (34) 
by means of which the terms of any sequence satisfying 
relation (33) can be expressed by formulas (39) (i. e. by 
multiplying the terms by certain numbers A, B, .... C , 
and adding the products) is called the basis of the recur¬ 
sion sequence. 

It follows from the above sequence that every recursion 
sequence has a basis which can be selected in different ways. 

For example, systems whose first terms are 

1,0, . . ., O'! 

0, 1, .... ol 

0,0, ..., 1, 

~ (k) 

form the basis of an arbitrary recursion sequence of the 
&-th order. 

Let us sum up the results of 5. 

For every recursion relation of order k there exists an 
infinit 1 number of different sequences which satisfy this 
relation. Any such sequence can be formed from k sequen¬ 
ces which satisfy the relation and form its basis by multip¬ 
lying each of the k-th sequences by certain corresponding 
numbers A, B, . . ., C respectively and adding the products 
termwise. 

Therefore in order to obtain a full solution of a recur¬ 
sion relation of order k it is sufficient to find but a finite 
number k of the sequences that satisfy it and form its basis. 

We now give some examples in order to make this for¬ 
mulation clear. 
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or 


1, 1, .... f 
0 , 1 , ..., 1 

0, 0, , 1 . 
(k) 



Example I. Let we have the following recursion relation 
of the second order: 

u„ +2 2«, i+1 u n . 

Its basis must consist of two sequences: 

x lt x 2 , x a , . . . , X n , . • • , 

Ui* yai • • ■ j Un* ' • • 

In selecting them we shall assume 

^ = 1, x 2 =l and i/i = 0, y 2 =\. 

Since rewriting the recursion relation in the form 

+ 2 ^n + l = U n +1 


we see that the difference of two successive terms of the 
sequence is constant, i. e. the sequence which satisfies this 
relation, is necessarily an arithmetic progression, then the 
sequence { x lt \ whose first terms are x 1 = 1 and x 2 = l is an 
arithmetic progression, with zero common difference, i. e. 

1, 1, 1, ..., 1, ... (x„=l) 

and the sequence {y n \ with y l — 0 and y 2 — 1 is an arith¬ 
metic progression with the common difference equal to 
unity, i. e. 

0, 1, 2. n— 1... (*/„ = n — 1). 

A term of any recursion sequence which satisfies this 
equation can be expressed in the form 

u„ = Ax n + By n = A+ B {n — 1) 

where A and B can be determined from the equations 
u x = A + 5(1-1), 
u 2 = A + B( 2 - 1 ), 
or 

u t = A, 

u 2 = A B- 

Hence 


A—u lt B — u t — Uj 

it follows that 


«„ = «! + (« — 1)(«2 —« i )- 
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This is the general formula for a term of any recursion 
sequence satisfying the relation 

U M + 2 “ 2u n + 1 U u . 

Denoting u 1 — a, u 2 — u 1 ==d we can give it the following 
form: 

u n = a-f (n— 1) d. 

This is the well known formula for the general term of 
the arithmetic progression. 

Example 2. Consider another recursion relation of the 
second order: 

U ,l + 2 = U ,l + 1 + ' 

Assuming x 1 = \, x 2 =l we obtain the familiar Fibonacci 
sequence: 

1, 1, 2,3, 5, 8, ... 

As a second sequence entering into the basis let us take 
the sequence {y n \ in which y l ~ 0 and y 2 = 1. We can write 
then 

f/3 = f/ 2 + t/i = l> t/4 = t/3+t/2 = 2, y, = y t + y, = 3, .. . 

In this sequence y 2 = x lt y 3 — x 2 , r/ 4 = x 3 , y 3 — x t .. . and in 
general y„—x n _ 1 (ti — 2, 3, .. .). Indeed, since we have found 
that these equalities hold for any value of so 

that in particular y m+1 =x m , y m = x m . 1 , then for y m+2 we 
shall have 

IJm + i dm +1 "T yrn ~ ~t~ X m — I — X m + j 

and the equalities also hold for n = m + 2. 

Thus 

Vn — -1 (ti — 2,3, ...). 

Therefore for any sequence, which satisfies the relation 
+ 2 = U/1 + 1 "t* , 

we find [formula (39)]: 

u n = Ax„ + By n , 

where A and B are found from the equations: 

Uj = Ax 1 -}~ By l = A , 
u 2 — Ax 2 By 2 = A -}“ B . 
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Hence 


and 


A = u 1} B = u 2 —«! 


u„ = u 1 x a + (u t — u 1 )y„. 

With n^2 x can be substituted for y n \ then 
u n = u l x n + (u 2 — u l )x n _ 1 (n>2); 


or 

“» = «l — + «.*„-1- 

With n>3 

*/i ~ *n -1 "F * ;i - 21 X,, *„-i = */i-2, 

and consequently 

u„ = w x x n _ 2 -f-u 2 -V„ _ t (n ^ 3). 

Thus terms of any sequence {«„}, which satisfy the re¬ 
lation 


+ 2 + 1 ~t~ W n , 

can be expressed in terms of Fibonacci numbers by the 
formula we have established. In particular if u 1 = — 3, 
« 2 =1 (cf. p. 16) then 

“ n = — 3 x „-. + (n>3). 


6. It will be shown now that under some very general 
conditions the basis of the recursion relation (33) 

^■n + k~ ^l^n+k-1 "F ^2^n + k-2 • • • *F &k^n’ 

comprising k geometric progressions with different common 
ratios can be found. In order to do it let us determine 
under what conditions the geometric progression 

*i = l, x % = q, . . ., x n = q n ~\ . . . (q¥=0) 

satisfies ralation (33). Noting that 

y - nil + A — 1 y - + 2 y .— nl 1“ 1 

x n + k — y > x n + k-\— H » •*•» x n~~H » 

we substitute these quantities in equation (33) (for u n+k , 
u n+k _ l , .... «„) and obtain: 

q n+k - 1 =ra 1 q n+k ~ 2 + a 2 q n+k ~ 3 + ■ ■ ■ + a„q n ~\ 
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whence 


q k = a x q k ~ 1 -f a % q k ~ 2 + . . . + a k . (43) 

Thus a geometric progression can satisfy the recursion rela¬ 
tion (33) of order k if and only if the common ratio q of 
the progression satisfies the algebraic equation (43) of de¬ 
gree k with the same coefficients as those in equation (33). 

Equation (43) is called the characteristic equation for 
the recursion relation (33). If q = a is a root of the cha¬ 
racteristic equation (real or imaginary) then assuming 

x n — cc n ~ 1 (n= 1,2,...), (44) 

we obtain a geometric progression whose first term is 
x x — 1 and with a common ratio a; this progression satis¬ 
fies relation (33). Since a is a root of equation (43) we 
have, indeed, 

a*== ajcc* -1 -\-ayx h ~ 2 -j- . . . +a k . 

Multiplying the two members of this expression by a" -1 
where n is an arbitrary natural number we obtain which 

a n+ ' ! - 1 a l a n+k ~ i + a 2 a n+k ~ 3 + . . . +a k a n ~ l , 

means that sequence (44) satisfies relation (33). 

Thus to each root q — a of the characteristic equation (43) 
there is a corresponding geometric progression (44) with a 
common ratio a; this progression satisfies the recursion 
relation (33). 

To form the basis of geometric progressions with 
different common ratios only one must have a sufficient 
number k of them and consequently have k different roots 
of the characteristic equation. 

Assume that all the roots of the characteristic equation 
are different: 

q x —(x, g 2 = p, . . ., <7& = Y> 

We have then k geometric progressions which satisfy 
relation (33): 


1, a, a 2 , . 


\ 

• » 



1, P. P 2 , 

P"-\ •• 

• » 

' ' 1 

> 

(45) 

i, v, y\ ■ ■ 

• , Y n "\ ■ ■ 

■ j 
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We shall show now that the system of sequences (45) is 
the basis of relation (33), i. e. that for any sequence {u n ) 
which satisfies relation (33) we can find numbers A, B,..., C 
such that for any n we shall have 

u„ = Aa"- 1 A-B$ n ~ x + ■ ■ • +Cy" -1 - (46) 


In order to prove this it suffices to verify that the sys¬ 
tem of equations obtained from (45) with n= 1, 2, . . k 


A-\- B • • • + C = u } , 

Act -|- Bp + • • • + Cy = u 2 , 


Aa k ~ l + B$ k -'+ . . . u k , ) 


(47) 


has a solution with respect to the unknowns A, B, . . C 
with any values of the right-hand members of these equa¬ 
tions; for this it is sufficient, in turn, (cf. p. 19) that 
the corresponding homogeneous system 

A + B+... +C = 0, 

Act + Bp + . . . -f Cy — 0, 


Aot k ~ x A- Bp A_1 -f . . . +Cy*- 1 = 0 



have one and only one zero solution. That is actually so. 

In fact, suppose there exists a zero solution of (48), i. e. 
there are numbers A, B, . .., C such that at least one of 
them, say, A is different from zero, and the numbers sa¬ 
tisfy system (48). To reduce this proposition to a contra¬ 
diction let us first form a polynomial M (x) of degreed—1 
which vanishes with x = p, ...,x — y and equals unity 
with x — a. Since this polynomial is of degree k — 1 and 
vanishes with k—\ different values of x: p, . . ., y, it must 
be of the following form: 

M (x) = \i(x— ^p) ... (-K--Y). 

where p, is some number. Taking x — a, we must obtain 
M (a) = 1; hence 


and 


1 = p (a P) . . . (a— y) 


_I_ 

— P) •••fa—Y)' 
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Thus 


M (x) = 


(* — P) (x—Y) . 
(« — P) • ■ ■ («—Y) ’ 


and obviously the polynomial satisfies the given conditions. 
Removing the parentheses and collecting like terms we can 
give the polynomial the following form: 

M (x) = m 0 + WjX + ...-)- m k _ 1 x k ~ 1 . 

If we multiply now equations (48) by m 0 , m,, .... m,,_ 1 
and add them term by term we obtain 

A (m 0 + /«!«+ .. . + tn k _ 1 a k ~ 1 ) + 

+ B {rn 0 -j- -f- •. . + +••• + 

+ C(m 0 + m 1 Y+ . . . + m k-iy k ~ 1 ) = 0. 
or 

AM(a) + BM($)+ . ■ ■ +CM (y) = 0. 

Since M (a)= 1, M(P) — 0, .... M (y) = 0 it follows that 

4 = 0, 

which contradicts our assumption. 

Thus system (48) has one and only one zero solution 
and system (47) must have a solution (one only) with any 
u 1 , u 2 , .. ., u k , and this, in turn, means, that system (45) 
is the basis of recursion relation (33). 

Thus we have found that for every recursion relation 

u n+k~ a i u n+k-lA- ■ ■ • +a k ll n , 

whose corresponding characteristic equation 

q k = a 1 q k ~ 1 + a 2 q k -*+ ■ ■ ■ +a k 

has different roots: q — a, q — §, ...,q — y, there exists a 
basis consisting of k geometric progressions with the com¬ 
mon ratios: a, p, ..., y. In other words for terms of any 
sequence {«„} which satisfies relation (33) there are k num¬ 
bers: A, B, ...,C [they are found from equations (47)] 
such that 

w„ = 4a"- 1 + i5p"- 1 -f .. . -f-Cy"- 1 (n= 1, 2, 3, ...). 
Let us sum up the ideas of 6. 
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A recursion relation of order k has a corresponding alge¬ 
braic equation of degree k with the same coefficients—its 
characteristic equation. Each root of the characteristic equa¬ 
tion is the common ratio of a geometric progression which 
satisfies the recursion relation. If all roots of the charac¬ 
teristic equation are different, k geometric progressions are 
obtained which constitute the basis of the recursion relation. 
It follows that in this case terms of any sequence which 
satisfies the recursion relation can be obtained by adding 
term by term certain geometric progressions (k in number). 

7. We now turn to the application of the results ob¬ 
tained. Consider first the Fibonacci sequence. Its recursion 
relation is as follows: 

+ 2 = U n + 1 "T u n 

and the characteristic equation (43) accordingly takes the 
form 

<7 2 = <7+ 1- 

Solving this equation we obtain two different real roots: 
a = j + j]/5 and P = yJ-'X 

The general term of the Fibonacci sequence can there¬ 
fore be written in the following form: 

u n = Aa n ~ 1 -f- Bf)"' 1 . 

In order to find the unknown coefficients A and B 
assume n— 1 and n = 2; we obtain 


u l =\ = A + B 

u 2 =\ = Aa + B$ = ±.(A + B) + ^(A-B). 


Solving this system we find: 


Hence 


A 


^5+1 

2 V'5 ’ 


£ = 


y~5 —i 
2 i-CT 


u 


n 


V 5+J 
2 y~b 


(\ + y~5 y- 1 , ii / ]_ |^5 y 

V 2 ) + 2 1/5 v 2 ) 


1 
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This is precisely the general expression for the Fibonacci 
numbers. At first sight the formula obtained looks unwieldy 
and poorly suited for calculations. However using it one 
can obtain many interesting results. We shall show, for 
instance, that the sum of squares of two adjacent Fibonacci 
numbers is a Fibonacci number as well. 

Indeed 



Hence 


«4+i + Wn = 



It follows that 
For example 


5+/5 , / H-/5 Y" 5-/5 1 
2 2 / 2 ) 

1+/5V" + 1 /l— /5Y" + H 

2 ) \ 2 ) f — U 2» + 1- 


u n+l + w n = U 2n+l (50) 


«i3 = “7 + «l= 13 2 + 8 2 = 233. 


By the way, this is the answer to the Fibonacci problem. 
We leave it for the reader to prove a more general rela¬ 
tion for the Fibonacci numbers than (50), viz. 

^■n^m “F + \U m + 1 = ^n + m + l- (^1) 


The following theorem will be proved as another appli¬ 
cation of formula (49): 

For any two natural numbers a and b with a < ft, the 
number of consecutive processes in the Euclidean algorithm 
which are necessary to find the greatest common divisor 
(G.C.D.) of a and b is not higher than five times the num¬ 
ber of digits in a written down using the decimal system. 
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Applying the Euclidean algorithm to find the G.C.D. 
of the numbers b and a we obtain a chain of equalities: 

(1) b = ax'-\-y', 

(2) a = y'x" + y", 

(3) y' = y"x'"+ y'", 

(52) 

(k) y (k ~ 2) = 

(fe-f 1) yC*-l> = ^*>^(*+l>. 

The consecutive remainders satisfy the inequalities 
a>y'>y"> y'" > . .. > z/'*- 1 ' > 1. 

The remainder in the last of equalities (52) is zero. 

This means that the foregoing remainder y M is the re¬ 
quired G.C.D. of the numbers b and a. Therefore k is the 
number of operations required for finding the G.C.D. Our 
problem, as we have said, consists in evaluating the num¬ 
ber k. For this purpose let us compare the numbers y {k \ 
y {k ~ 1} , .. ., y', a with the Fibonacci numbers u lt u 2 , u 3 .. . . 
Note that y {k) ^\ = u 2 but the foregoing remainder 
is greater than y ik) , and hence # (A-1) ^ 2 = « 3 . Therefore 
we conclude from the equality k that 

y {h ~ 21 = y (R ~ l) x {k) + y {k) ^ y ik ~ 11 • 1 + y u<) ^ » 3 + = « 4 

Thus 

if k) ^u 2 , y ik ~ 1] ^ « :i , y {k ~ 2) ^u i . 

Assume that we have proved the following inequalities 
to be valid 

y {l<) >u. 2 , ..., y im) >u k _ m+2 , y {m ~ 1) >u ll _,„_, 3 
(m—1 >2). 


Then from the equality y Un ~ 2) = y {m ~ l) x m y m we obtain 

y un- =) ^ (/<»-!>. 1 + y^n) ^ Uk _ m+3 _j_ Uk _ m + 2 = Uh _ m+4 . 

Continuing this argument we come to the inequalities 

y" y’>Hk +1 

and further from (52, 2) we deduce that 

a—y'x" + y" ^ y' ■ 1 + y ^ w *+i + u / ; — u *+2- 
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By formula (49) the term u h+2 takes the form 


« + 2 " 


V 5 


14-/5 


<r \*+2 


-(^H 


Therefore 


a>- 7 = 

/5 


1 + /5 


ft +2 


1- /5 V 


> 


(since 


V6 


> / 5 


i r (i + /s 


r\*+! 


— 1 


( 53 ) 


< 1, and consequently 


1 — /5 


k + 2 


<1 • 


From (53) we find that 

l±V±y‘ <aVl+ 1 < K5(« +1) < (!+/!)■(«+1) 


r, /i+rjY = 3 + y~5 


(K5<( 


Therefore 

(■ 

We now note that 

“• = 6 “ 7 ?[ 

Hence 
Consequently 

10* < 


J 


+ /5 


as 


/5<3) 


•< rz 1. 


(54) 


17 1+/5 V 


IA 2 ) 

V 2 )\ 


< 


< /5 


1+ /5 \ 5 


+ 1 


1 + /5 


>5/5—1 > 10. 


1+ /5 


5ft 


</+!)> 


(55) 


If number a is written down in the decimal system with n 
digits (a is an n-digit number) it is obvious that 

10" _1 < 10"; 

hence 

a+ 1 < 10" 

and thus by virtue of inequality (55) 

10*<(a+l) 6 <10 5 " 
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or 

« < 5 n. (56 

This is precisely the required result: the number of con¬ 
secutive k divisions in the Euclidean algorithm is less than 
five times the number of digits in the smaller of the two 
numbers, b, and a, written in the decimal system. The 
foregoing proof shows that the most unfavourable case of 
the application of the Euclidean algorithm (as to the num¬ 
ber of operations near the limit established by the theorem) 
is that in which b and a are the adjacent Fibonacci num¬ 
bers. In order to confirm this statement let us take, for 
instance, b — u. 20 = 6765 and a = Ui 9 = 4181. Here a is a 
four-digit number and by the above theorem the number 
of operations in the Euclidean algorithm must be less than 
5x4 = 20. In fact we have to perform in this case k = \7 
operations: 

(1) 6765 = 4181x 1 + 2584 (10) 89 = 55x 1+34 

(2) 4181 = 2584x 1 + 1597 (11) 55 =34x 1+21 

(3) 2584= 1597x 1 + 987 (12) 34 = 21x1 + 13 

(4) 1597= 987x1+ 610 (13) 21 = 13x 1 + 8 

(5) 987= 610x1 + 377 (14) 13= 8x 1+ 5 

(6) 610= 377x1+ 233 (15) 8= 5x1+ 3 

(7) 377= 233x1 + 144 (16) 5= 3x1+ 2 

(8) 233= 144x1+ 89 (17) 3= 2x 1 + 1 

(9) 144= 89x1 + 55 (18) 2= 1x2+ 0 

Fibonacci numbers in descending order are obtained here 
one after another as the remainders of division. All quoti¬ 
ents (except the last) are unity and this is why the number 
of operations is so great. The greatest common divisor 
proved to be unity [equality (17)] which could be foreseen 
beforehand for the adjacent Fibonacci numbers. Indeed, 
from « rt+2 = + «„ it follows that the G.C.D. of the 

numbers u n+2 and u n + 1 is the same as the G.C.D. of u n + l 

and Therefore for any pair of adjacent Fibonacci num¬ 
bers the G.C.D. is the same. In order to find it one can 

just consider the pair u 2 = u 1 -= 1 which shows that the 

G.C.D. is unity. 

8. Let us consider as the next example the recurrent 
sequence (16): 

«! = 5, u 2 = 7, u 3 = 1, w 4 = 3, m 5 = 2, u 6 =l, «, = 3, ... 
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The recursion relation in this case is 
««+3 = (n>3), 

and thus the characteristic equation is 

q 3 = 1. 

The roots of this equation are 


.. ITJ 


a=l, p = — + antl Y — — t —* 


2 1 2 » — 2 

Therefore the general term of the sequence should be as¬ 
sumed to have the following form [cf. formula (46)]: 

u n = Aa"- 1 + Bp" -1 + Cy n ~ 1 = 
l , . / I '"' 1 


— A -f- B ( —g- + i 


+ C 


lYLY -1 

2 2 y 


We can require this formula to be valid for all values 
of n with which our recursion relation holds too: n = 3, 
4, 5, .... 

Note that 


JL + i^± 
2 ‘ 1 2 


1 


2 1 


, Y3- 


JX > . 3X 

COS y— l Sin y 

31 . . . 3X 

COSy + J Siny 


+ 


Therefore by De Moivre’s theorem (formula) 

= /4 -J-(—1)" -1 B cosy (n— 1) — i siny (n — 1) 

+ (—l) n_1 C |cosy (n — 1) + J siny («— 1) = 

= yi + (B + C)(—l)"- 1 cosy(rt—l) + 

+ t -(-5 + C)(-l)'’- J sinf (n-1). 

Let us have B-\-C = A v and i (— B + C) = A 2 ; then the 
above formula will be rewritten in the form 

u„ = A + (— l)" -1 cosy (n — 1) + A 2 (— l)"- 1 siny (n — 1) 

(«> 3), 
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and we have now to find only the unknown coefficients 
A, A x and A 2 . Taking n — 3, n — 4 and n — 5 we obtain 
three equations in three unknowns: 

A-^A. + ^-A.,, 

A + A lt 

-A-^A.-^-A,. 

Hence we find 

^4 = 2, A 1 = 1 and A 2 = — -y=- . 

Thus 

«» = 2 + ( —l)”- 1 COS(/l— 1)-J — yLrsin(« — l)y = 
^2 + (-l) n ~sin(n-2)f (»> 3). 

It can be seen that in this example the general term of 
the sequence is expressed by trigonometric functions which 
is in agreement with the recurrent character of the sequence. 

As the last example consider the one directly related to 
the division of polynomials. 

Given two polynomials: P (x) — 3 + x 2 —x 5 and Q (x) = 
— 2—x—2x 2 -f-x 3 . The problem is to find the structure of 
the coefficients in the quotient of P(x) by Q (x). The se¬ 
quence of the coefficients in the quotient 

u l — D n , u 2 — D l , ..., u n — D n _ j, ... 

as noted in 2, is a recursion sequence whose terms satisfy 
equation (20): 

D n+k = - |f D n+k _ x — ... -|f D n (n > l-k + 1), 

where k is the degree of Q (x), B 0 , B x , ...,B k are the 
coefficients of Q (x) and Z is the degree of P ( x ). 


u 3 — 1 = A + A x cos -|- A sin = 
u 4 = 3 — A — A x cos^-— A 2 sin ^ - 
u 5 = 2 — A +A X cos -f A 2 sin y- = 
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This means that k — 3, B 0 = 2, B 1 = —1, B 2 — ~ 2, B a = 1, 
1 = 5 : 

Dn+a — y Ai+a + y A>+i y B>„ (n^5 3+1=3), 

i. e. 

B^n + 3 ~ y D n+2 + D n +i ( n ==^^). 

The characteristic equation takes the form 

? 8 = y?* + ?~ y> 

or 

q*-q-±(q*- 1) = (q-\) (q- 1) (q + 1) = 0. 

Therefore its roots are 

a = y > P= 1. Y = —1. 

and for D n we obtain the formula 

D n =^(4) n + fl.l» + C(-l)" (n> 3). 

Assuming n = 3, n = 4, and n = 5 we obtain three equa¬ 
tions: 

D a = -g- A -j- B — C, 

Dt — jftA-j-B-^C, 

d & =±a+b-c. 


In these equations not only the coefficients A, B, and C, 
but also the numbers D 3 , D t , Z) 6 are unknown. In order to 
find the unknowns let us perform the division of P (x) by 
Q (*) up to the point where the quotient will contain terms 
of the fifth power, inclusive. 
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Dividing we get 

3-f-A 2 — x s 12—a—2a 2 + a 3 _ 

o ^ 0 „ . 3 „ 3.3 , 0 3 „ ,3 , r.19 . 

3 yA 3 a + yA ~2~\~T X ~^ Ti* + 

yA + 4A 2 —yA 3 —A 5 
3 3 „ 3 , . 3 . 

~2 X ~\ x ~2 X ~^T X 

. 3 , 3 . . 

_ 4 4 

4-| A 2 —- 2y A 3 — 4y A 4 2-|- A 6 

_4* i + 4 * , -4*‘ 

2 it jc;i ~*T 6 j:4 ~ 2 '8 * b + *T6* 6 

r- 3 . . , 3 „ 

5 s*‘~ 2 i*‘- 5 A*'+ 2 l*’ 

1$ 


Hence 




D = 2— 

U l L s > 


D. = >ra 




and the system of equations obtained above takes the form 

i 4 + B + c - 2 I' 

whence we find 

A—4— B-— C- — 

Thus 

= 4y • gs + T + 1T ( rt ^3). 

The problem is solved. It follows by the above formula 


D = 2— D = — 

^ 1 98 > ^7 ocfi ' 


£) — 2 — 

8 z 512 ’ 
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9. In all of the above examples the characteristic equa¬ 
tion had only simple roots. Now consider the example 
pertaining to the sequence whose terms are the sums of 
squares of natural numbers (given on p. 15). The recursion 
relation for this sequence is 

Ni +4 = 4s „ +3 6 s n + 2 + 4s „ +1 s„, 

and accordingly its characteristic equation takes the form 
q 4 = 4 q 3 — 6 q* + 4q — 1, 
or 

q* — 4 ( 7 3 + 6 < 7 3 — 4 q-\- 1 = (< 7 — l ) 4 = 0 . 

It has only one quadruple root: < 7 = 1 ; therefore we ob¬ 
tain in this case only one geometric progression with the 
common ratio 1 whose terms satisfy the above recursion 
relation. 

In such cases we have to find other simple recursion 
sequences which can form, together with this geometric 
progression, the basis of our equation. Such sequences in 
our example are: 

0, 1, 2, 3, 1, ...; 

0, 1, 4, 9, .... (<i-l) 2 , ...; 

0, 1 , 8 , 27.(»-1)». 

(the reader can easily verify this assertion). Without treating 
this very general case which would require rather unwieldy 
calculations we shall consider the following typical example. 
Let the recursion relation be 

Un+k = Cfc Xa *n+k-i— c k- i a?u n+ k- 2 + ••• +(—1 ) k ~ l C a k a k u n , 

(57) 

where C| _I , Cl -2 , ...,C| are binomial coefficients of the 
6 -th order. The corresponding characteristic equation 

q k = C k k -'aq k -'— C|-*aV"*+ . . . +(—lJ^Cgo* 

can be written in the form 

(i q — a) k — 0. 

It has a 6 -fold root <7 = a; obviously 
(a—a)* = a*—Cfr^a* + C£~ 2 a*— .. . + (— \) k C° k a k = 0. (58) 
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Consider in general the following identities 
(a— a) k ~' n = a k - n — Cfc'^ 1 a k ~ m + 

+ Cfc*-*a fe -«- .. . + (-1 f- m C%- m 0L k -* = 0, 

where m — 0, 1, 2, ..., £—1 or 

(1 - \) k ~ m = CtZ-CtZ- 1 + • • • + 

+ (-l)^:r+ • • • + (—l) ft_m CL,„ = 0. (59) 

Equality (59) corresponding to m = 0 takes the form 

C*-CJh + Cjb-. • • +(-!)•* C{-“ + 

+ ...+(-i) fc a=o. ( } 

Noting that 

l) ... (p + l) fe(fe— l) ... (*— m+\) rk _ m _u 

_ 1-2 ... (fe — n) ~(k — m — p+1) ... (* — p) k ~ m 
(m— 1, 2, ...,k —1; —m), 

or 

k(k-l)...(k-m+l)Cil*-» = 

==(k — m-—\i,-\-\)...{k — \i)C k k~ v ', (60) 

we multiply each of the equalities (59') by the correspon¬ 
ding multiplier k(k —1) ... ( k —m+1). We can now using (60) 
write them in the form 

(. k—m + 1) ... kCk — (k — m) ...{k— 1 )Cjb+ ...+ 

+ (—1) M '(^— m —(* + 1) • • • (k — 

+ ... -f (—l) fe - m l-2 ... ma = 0 (59") 
(m= 1,2,..., k —1). 

We shall now prove that the following equations hold 
with m = 0, 1, 2, ..., k —1: 

k*ci- (k -1)»a- 1 +...+(-1) •* (k-nr ct »+ 

+ • • • + (—l) ft -0'”Cft = 0. (61) 

Indeed, the equality corresponding to m = 0 coincide 
with (59') and thus holds true. 

Assume, reasoning by induction, that equalities (61) have 

already been proved for m —0, 1./ (j^k —2). Let 

us demonstrate that the equality which corresponds to 
m = j +1 also holds true. In order to prove it we shall 
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make use of a polynomial of degree /-fl: 

f(x) = (x—j) (x—j+l).. .{x— l)x = 

= xi+'—fyx !—.. . — fax. (62) 

Multiplying (61) with m— 1,2, ..., j by the correspond¬ 
ing numbers (3^ P 2 , ..., fly we obtain 

p 1 *ci-M*-i)c*- i + 

+ ...+M-D* (A-ji) c k k ^ + 

+ • • • + Pi (—l)*-0-Cft = 0, 

. (63) 

V J vct-b(k-iyci- i + 

+... +p y (—l)**- (k-iiyc k k -»+ 

+ • • • + Py(—l) ft -0-Cfe = 0. 

Let us rewrite equality (59") corresponding to the value 
m = / + 1 in the form 

/ (k) a-f (k - 1) a- 1 +..,+(-1 rnk-v) + 

+ ...+(-l)*/(0)a = 0. (64) 
In doing that we made use of the equalities 
(k-j)...k = f(k), 

(k-l-l)...(k-l) = f(k-l), .... (k-VL-j)... 

■ • - (k — j) = f{k —/),— 

Adding (63) and (64) term by term we have: 

[Pi&+ • • • +P jk J J rf {k)\ C\ — 

—[Pi a* —i) +... +P j{k —i y ~t~f (k —i)] + • • • + 

+ (-1)^ [Pi (As—I*) + • • • +Py (k-iiY + Hk-ii)] c$r» + 

+ ... +(—l)*[Pi-0+ ... +p,.•()'•+/(0)] C% = 0 

But by (62) we can write 

Pi* 4- Pa* 2 + • • • + Py* y 4* f (*) = x J+1 . 

Therefore the above result yields 

v+'tt-{k-\y+'cfr'+ ...+(-iy* {k- v ) i + 1 cl» + 

+... +(—i)*-o/ +1 -cx = o. 

This is precisely equality (61) for m = j + 1. Thus we 
have demonstrated that Eqs. (61) hold true. 
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Lastly, consider an arbitrary polynomial whose degree 
is not higher than k — 1: 

P (x) = A k _,x*~' + A k _,x^ + ... + A,. (65) 

Multiplying equalities (61) with m = 0, 1, 2, fe — 1 
by corresponding factors A 0 , A lt .... A k _ 1 we obtain: 

A 0 C k — A 0 Ck~ l -\- ••• + 

+ (-l)M 0 C^+...+(-l)M 0 CJ = 0, 
A 1 kC k k -A 1 (k-\)C k k - l +...+ 

+ (—\Y A 1 (k — p.) C£ _M, + • • ■ +(—1)* ^i'0-Cft = 0, 


A k _^~'Cl-A k _, (k-l ...+ 

+ (-l)' i A k _ 1 (k-\i) k ~ l C k k -»+ ... + 

+ H)M*. 1 -0‘- 1 .C| = 0. 

Adding the equalities term by term we obtain 

(^o+ Ajt +...+ A k _ 1 k k ~ i ) a — 

-[A 0 + A 1 (k-l)+.--+A i _ 1 (k-l)*~ 1 } C|- 1 +...+ 

+ (-\r [A 0 + A, (k-ii) + ... + 4*-i (^-P) 4 - 1 ] Ci-“ + 
+ ...+(-l)n^„ + A 1 -0+...+A A _ 1 .0^ 1 ]Q = 0, 


or 


P (*) • Cl- P (k -1) ■ cl 1 + ... + (-1 ) k P (0). ci = 0. ( 66 ) 

And hence an arbitrary polynomial P (x) whose degree 
is not higher than k —1 satisfies relation (66). 

Assume, in particular, P(x) = (x + n —l) m , where n is 
an arbitrary natural number and m, an integer 
— 1 . 

Then equality (66) takes the form 

(k + n-l) m C k k -(k + n-2rC k k - 1 +. .. + 

+ (-l)*(n-l)"C* = 0, 

or multiplying it by a k+n ~ 1 and substituting unity for C*: 

(H«-l)'"a H "' 1 = C|' 1 a(H«-2)" , a' i+ "- 2 - 
—C k k ~ 2 a 2 (k + n— 3)”a* +n ~ 3 + ...+ 

+ (— l)*- 1 Cga* (»—!)" a"- 1 . (67) 
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Comparing (67) with (57) we see that each of the fol¬ 
lowing k sequences satisfies the recursion relation (57): 

1, a*, ....a"' 1 . (m = 0); ^ 

0, a, 2a 2 , (n —l)a" _1 , .... (m= 1); 

0, a, 2 2 a 2 , (n —l) 2 a n_1 , (m = 2); > (68) 

0, a, 2*~ 1 a 2 , .... (n — l)* -1 a 0-1 , .... ( m — k — 1). j 

If we come to establish that these equalities form a ba¬ 
sis it will follow that the general term of any sequence 
which satisfies equation (57) has the form 

u n = [B 0 + B 1 (n—\) + .. .+B k _ 1 (n—\) k ~ 1 ] a"' 1 = 

= Q (ft — 1) a" -1 , (69) 

where Q (x) — B 0 + B 1 x+ . .. J r B k _ 1 x k ~ 1 is a polynomial of 
a degree not higher than k — 1 with arbitrary coefficients. 

In order to prove that (68) is a basis it is sufficient to 
demonstrate that the following system of k linear equations 

B 0 4* Bi ■ 0 4~ • • • 4" B k _ x • 0 A 1 = u 1 , 

5 0 + B x - 1+ ... -\-B k _ 1 -1 4-1 = w 2 > 


B 0 + B 1 (k-l)+ ■. ■ +B k _ 1 (k-l) k ~ 1 = u k 

has a solution with respect to the unknowns B 0< B lt .... B k _ 1 
with any n 1( ...,«* (which follows from the proposition 
stated on p. 19), i.e. that the system 

B 0 — 0, 

B a + B x -f-... -f- B k _ i — 0, 

B 0 + {k-l)B l +...+(k-l) l ‘-'B k _ l = 0 

has only a zero solution. But the equations of the latter 
system show that 

Q(0) = Q(l)=-..=Q(*—1)==0, 

i.e. the equation 

B 0 -\-Bj^x-}-B^^ x x k 1 = 0 

of a degree not higher than k —1 has at least k different 
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roots: 0, 1, 2, .... k —1. Hence 

B Q = B 1 = . .. = B k _ 1 — 0, 

and thus we have proved that sequences (68) form the 
basis of the recursion relations which satisfy equation (57). 
If the recursion sequence satisfies the general equation 

“n + k = a iU„ +k - 1 + a 2 u n+k _. 2 + . . . + a A u n (a, { =£ 0), (70) 
its characteristic equation 

q k = a v q k ~ 1 + • • • + a k (71) 

can have a root a of multiplicity a, a root p of multip¬ 
licity b, .... a root y of multiplicity c, and a total of 
a + b + ... +c = k roots. 

In this most general case it can be demonstrated that 
the basis is formed by the following k sequences: 

1, a, a 2 , .. ., a"" 1 . 


0, a, 2 a ~'a 2 , .... (n — l)“ _1 a' ,_1 , ...; 
1, p, p*. ..., p"~\ .... 


0, p, 2 6-i P 2 .(/i-1) 6 " 1 P"- 1 , ...; 


1. V. Y*» • • •• Y" 1 » • • •» 


0, y, 2 c-1 y 2 , .... (n— l) c_1 Y n_1 , ... 

Therefore 

«>*— Q{n— Ija"- 1 R(n— l)P n_1 + ... +5(« — l)v" _1 , (72) 

where Q(x), R(x), ..., S (x) are some fixed polynomials 
of a degree not higher than a—1,6—1, ...,c —1, respec¬ 
tively. 

It follows that the general term u n of any recursion 
sequence has the form of a sum of the products of polyno¬ 
mials with respect to (n—1) (or, which is the same with 
respect to n) by general terms of geometric progressions 
whose common ratios are equal to the roots of the characte¬ 
ristic equation (71). 

If all the roots of the latter equation are simple, the 
above polynomials are constants and the general term of 
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the recursion sequence takes the form of a sum of the terms 
of geometric progressions. 

It can be proved that the converse also holds true, na¬ 
mely, that any sequence \u n \ whose general term is express¬ 
ed by formula (72) is a recursion sequence*. 

The corresponding equation (71) is constructed using its 
roots a, p, .... y and their multiplicities a, b, ..., c (which 
are the increased by unity degrees of the polynomials Q, 
R, Hence the recursion relation (Eq. 70) can be 

found immediately. 

Consider as an example the sequence 
u„ = (n — l) 2 .2 n-1 + 3 n-1 

Comparing it with (72) we see that the roots of tne 
characteristic equation are: a = 2, (3 = 3, a being a mul¬ 
tiple root of multiplicity 2+1 = 3. Therefore the charac¬ 
teristic equation must take the following form: 

(q— 2) 3 (q— 3) = q* — 9q* + 30? 2 — Uq + 24 = 0, 
and the recursion relation is 

u n+i — 9 m„ + 3 — 30« n+2 + 44«„ +1 — 24«„. 

It is left for the reader to ascertain that the above se¬ 
quence satisfies the last equation. 

10. We shall illustrate now the results of section 9 by 
several examples. It was shown in section 2 that the terms 
of an arithmetic progression satisfy the relation 

^n +2 2W„ +1 U n , 

the squares of natural numbers, the relation 
^n + 3 = 3u, i + 2 3 u n + l + u n , 

and the cubes of natural numbers, the relation 

u n + i = 4w n + 3 + 2 + 4 U ll+1 U n . 

Obviously all these equations are particular cases of the 
equation 


* The proof of the theorem and its converse, which we omit, can 
be found in Chapter 4 of the book Division which Leaves a Remainder 
in Arithmetic and Algebra, Moscow-Leningrad, 1949; the theory of 
recursion sequences is expounded in this book in a manner different 
from that used in this booklet. 
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u n+k = Ct 1 u n+k - 1 -C k k -*u n+k _ 2 + ... -M-l)*- 1 C|«„, (57') 

discussed in section 9 (in this case a=l). 

The general term of any sequence which satisfies this 
equation must have the following form [cf. formula (69)]: 

u„feS 0 + 5 1 (n-l) + ...+5,_ 1 (n-l)*- 1 . (69') 

In order to find the coefficients B 0 , B lt ... B k _ 1 it is 
sufficient to solve the following system of k linear algebraic 
equations in k unknowns: 

B 0 = Wj , 

B 0 + B 1 + • • • + B k _ 1 = u 2 , 


B 0 + B, (k — 1) -f- ... + B k _ 1 (k — 1) A 1 — u k . 

In the case of an arithmetic progression 6 = 2 formula 
(69') takes the form 

= + B x (n —1), 

and system (73), the form 

B a ~u l , B n -\-By — u 2 . 

It can be seen that B 0 = u 1 is the first term of the pro¬ 
gression and B t = u 2 — u 1 = d, its common difference. It fol¬ 
lows that 

u n = u 1 + d(ti— 1 ) 

which is a well known formula. 

There is no need in making calculations for the case of 
sequences of squares or cubes of natural numbers since we 
know beforehand that the general terms are u„ — n 2 or 
u n = ri J . However, there is some interest in applying the 
relations (69') and (73) to the derivation of formulas for 
the sum of an arithmetic progression, as well as for the 
sum of squares or cubes of natural numbers. 

In section 4 we proved that if the terms of a sequence 
\u n ) satisfy the equation of the form 

u n+h ~ a i u n + k~l a 2 U n + k-2 4" • • • "t~ a k u m 

then the sums {s„} of the terms of these sequences (s^Mj, 
s 2 = u y -\- u 2 , s 3 = u t + u 2 + u 3 , ...) satisfy the equation 

• S n + * +1 = (1 + fl l) S n+k + ( a 2 - a i) S n + /{-1 + • * • + 

+ i a k — a k-l) S n+l a k s n- 
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Obviously for equation (57') we have 

= a 2 = -Cl ...,a k = (-l)*~'C k k . 


Therefore 

1 + a i = 1 4~ Cl = CL,» 

« 2 — (C| -|- CL -C| + i» 

q 3— a 2 = C* + CJ = CJ+i, 

a* - a* -1 = ( -1 )*" 1 (C| + C^) = (-1 )* “ - 


— a k = (—1 )* Cft = (— 1 )* Cft+i 

and the equation for {s„} can be written in the following 
form 


S n + A + 1— C\ +1 S, ) + t C* + ]S„ + ft -i 4* ■ • • ~T ( l)*C*tl S H> 


or 

S n+A +1 ££+l S n + A 4" CI + 1 S /! +A- 1 - ~• ■ • H~ ( — “1 )* + 1 Cfc+l S n = 0- 


Thus, if a sequence {u n \ satisfies the recursion relation 
(57') of the &-th order, then the sequence of the corres¬ 
ponding sums {s„} satisfies a relation of the same form 
but of order &+1- In particular, for the arithmetic pro¬ 
gression k = 2, for the sequence of squares of natural num¬ 
bers k — 3 and for the sequence of their cubes k = \ and 
it follows that for the sequences of the respective sums we 
must take k in the above equalities (57'), (69'), (73) great¬ 
er by unity: 3, 4 and 5. 

(a) The Sum of the Terms of an Arithmetic Progression. 
According to what has been said s n is expressed by for¬ 
mula (69') (substituting s„ for u n ) with k = 3. It follows 
that 

s n —B o + ^i ( n —1) + ^2 (n— 


The coefficients B 0 , B u B 2 , are dete rom system 

(73) (also substituting s„ for u n with k 


B o — Sj — u lt 

B 0 -f- B t -f- B 2 = s 2 = u x -(- u 2 = 2mj -f- d, 

B 0 -f- 2 • -f~ 2 2 B 2 — s 3 = Uj + « a + u 3 — -f 3d. 


Solving this system of equations we obtain 
B 0 — « 1( = r/j: -1 ^d, B% — d. 
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Whence it follows that 

s« = «i + («i l) + 4 d ( tt — 1 ) 9 = 

=««, + ( /z _ 1) n = 

_«[WH-«! + («— l)rf]_n («! + «„) 

2 ~ 2 

(b) T/ie Sum of Squares of Natural Numbers. 

Taking in formulas (69') and (73) & = 4 and substitut¬ 
ing s n for u n we have 

— ( n — 1) -\-B 2 (tx — l) 2 -f- B 3 (n — l) 3 

and 

B 0 = Si = 1» 

B 0 ~T B v -f- B. t -T B.j — s 2 — 1 -|- 2 2 = 5, 

B 0 ~T 25j -f- 4 B. t -f- 85 ;) = s a = 1 -f- 2 2 -f- 3 2 = 14, 

B 0 -f 3By + 9j3 2 27B 3 = s 4 = 1 + 2 2 -f- 3 2 + 4 2 — 30. 

The last system yields 

£ 0 =1, Bx = 2-i, 5 2 = 4> = 

Therefore 

s n =l + '-§(n-l) + ±(n-l)*+±(n-l)* = 

1 .. , 1 .., , 1 ..3 n(l+3« + 2n*) n (n+l)(2n + l) 

--6 rt + T n +T n --" 6-~-6-• 

We have obtained a well known formula. 

(c) The Sum of Cubes of Natural Numbers. This sum is 
expressed by the formula 

„ .«*(«+1) 2 
n~ 4 

We leave the derivation of this formula for the reader 
as an exercise. 

In conclusion let us consider one more example, the 
sequence a, 2a 2 , 3a 3 , . .., na n , . .. (a^O, a=^ 1). 

In this sequence 

u n — ria n (n — 1,2,3, ...). 

II is easily seen that 

— 2a«„ ,, —a 2 «„. 
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Indeed 


2 au„ + i—a 2 «„ = 2 a («+ l)a n+ 1 —a 2 raa n = 

= (n-\-2)a n+2 = u. n+2 . 

Since k = 2, a 1 = 2a and a 2 = — a 2 , the sequence consist¬ 
ing of sums {s„} (s 1 = a, s 2 = a-f- 2 a 2 , s 3 = a + 2 a 2 -f- 

-f 3a 3 , ...) satisfies the equation [cf. (30)]: 

s„ + 3 = («! + 1 ) s n+2 + (a. — a,) s n+1 a 2 s n = 

= ( 2 a + 1 ) s a + , — (a 2 + 2 a) s„ +x + a 2 s„. 

The corresponding characteristic equation is as follows: 
q 3 = ( 2 a-f 1 ) q 2 — (a 2 -f 2 a) <7 + a 2 . 

It can be easily ascertained that this equation is satisfied 
with q = a. Dividing the polynomial < 7 3 — (2a + 1) < 7 2 -|- 
-f- (a 2 4 - 2 a) q —a 2 by q — a we obtain the quotient 

q 2 — (a+ l)q + a. 

For this reason the two remaining roots of the characte¬ 
ristic equation satisfy the following equation: 

? J -(a+l)? + a = 0. 


These roots are a and 1. 

Thus the characteristic equation has a multiple root a 
of multiplicity a = 2 and a simple root {3 = 1. 

Therefore we obtain for s n [cf. formula (69) in which we 
must write s„ instead of u n , a = a, Q (x) = B 0 + B x x, i.e. 
a polynomial of degree 1, p = l and R(x) — C 0 , a constant]: 


s„ = [£o + £i (rt—l)]cc n ‘T-Co (n= 1,2,3,...)- 


The coefficients B 0 , B 1 and C 0 are found from the fol¬ 
lowing system of equations with n = l, 2 and 3, respecti¬ 
vely: 

B 0 + = s i — a > (Bo + Bj)a -f- C„ = s 2 — a -j- 2a 2 , 

(. B 0 + 2 Bj) a 2 -f- C 0 = s 3 = a -j- 2a 2 + 3a 3 . 

Hence 


B„ 


a 3 —2a 2 
(a— l) 2 


B x 


a 2 

a— 1 


and 


C 0 


a 


Consequently 

s n — [B 0 + B x (n l)]a." 1 -|-C 0 = 


na n + 2 — (n+l)a n + 1 +a 
(a— l ) 2 


u n a 2 — (u n+1 — u n ) 

(a-1) 2 
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CONCLUSION 


This book is meant to give the reader an idea of the 
variety of recursion sequences and of their place in mathe¬ 
matics. Along with it it was shown that the recursion 
sequences do not greatly differ from their simplest forms — 
the geometric progression and the sequences whose terms 
are powers of natural numbers (in particular, the sequences 
of natural numbers as such being an arithmetic progres¬ 
sion); the recursion sequences can be expressed making use 
of these simplest sequences. 

However even in elementary mathematics sequences are 
very often met with which are not recursion sequences. 
One of these is the sequence of prime numbers, one of the 
most important in mathematics: 

2, 3, 5, 7, 11, 13, 17, 19, 23, ... 

This sequence with its deep and complex properties is 
a subject of the theory of numbers. 

Sequences of values of many elementary functions 
are not recursion sequences as well. For instance, the 
sequence 

. _i_ j_ _i_ ± 

I, 2 ’ 3 ’ 4 ’ ’‘ n ’ '• 

(this is the sequence of values of the function y — ^ with 
x = l, 2, 3, ...), or the sequences 

i, V2, V3, y~4, .... vn, .... 

log 1, log2, log3, log4, .... logn, ... 

(these are the sequences of values of the functions ]/x and 
log jc) etc. 

These and similar sequences * (including the recursion 
sequences) are the subject of the branch of mathematics 
mentioned above — the calculus of finite differences. 


* Sequences of values of the so-called analytical functions are 
meant whose simplest instances are the elementary functions. 
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And finally, in elementary mathematics and particularly 
in analysis studied in colleges and universities a very 
important place is occupied by convergent sequences, i. e. 
sequences which have finite limits. Their study is the most 
important problem in the theory of limits and belongs to 
the foundations of analysis. The properties of individual 
terms of the sequences are even less than secondary: only 
the existence of a limit and its value are of consequence. 

We consider these remarks to be necessary in order to 
bring home to the reader that the theory of recursion 
sequences given in this book from the viewpoint of its 
subject-matter as well as from the viewpoint of the rela¬ 
tionships that result from it is only a particular and very 
modest chapter in the theory of sequences. 
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